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EXECUTIVE  SUMMARY 


This  two-part  paper  is  concerned  with  the  application  of  stereographic  projec- 
tion to  long-range  surveillance  systems  employing  a multitude  of  radars.  It 
is  an  extension  of  earlier  investigations  of  the  same  subject  in  connection 
with  the  Semi-Automatic  Ground  Environment  System  (SAGE)  and  the  National 
Airspace  System  (NAS).  Reconsideration  of  the  topic  at  this  time  is  appro- 
priate in  view  of  current  trends  in  the  air  traffic  control  (ATC)  community. 
Among  these  are  improved  accuracy  of  modern  sensors  now  being  contemplated  as 
replacements  for  present  operational  radars,  and  the  necessity  of  providing 
reliable  forecasts  of  target  position  over  long  time  intervals  in  advanced  ATC 
functions  such  as  Conflict  Alert.  A fundamental  consideration  is  whether 
errors  introduced  by  system  implementation  of  stereographic  projection  are 
compatible  with  the  quality  of  data  offered  by  new  radar  technology  and  the 
objectives  of  advanced  ATC  services.  Accordingly,  this  paper  deals  with  the 
origin  of  such  errors,  the  size  of  the  errors,  and  alternative  methods  of 
controlling  the  errors  in  the  context  of  the  reference  ellipsoid,  a mathe- 
matical representation  of  the  equipotential  surface  of  the  gravity  field  at 
mean  sea  level. 

In  long-range  surveillance  systems,  target  positions  are  represented  by 
points  in  a so-called  master  plane.  These  points  are  obtained  by  mapping 
target  altitude,  azimuth,  and  slant  range  into  the  master  plane.  If  the 
mapping  is  equivalent  to  stereographic  projection  of  target  latitude  and 
longitude  into  the  plane,  then  the  master  plane  representation  of  a target  is 
invariant  to  the  radar  from  which  azimuth  and  slant  range  are  measured  (pro- 
vided, of  course,  that  there  is  no  measurement  error).  In  a strict  theo- 
retical sense,  it  is  possible  to  achieve  such  an  equivalence  by  means  of  a 
two-stage  procedure.  First,  there  is  a conversion  of  target  altitude  together 
with  target  azimuth  and  slant  range  relative  to  any  given  radar  site  into  a 
point  on  a plane  unique  to  the  radar.  This  is  followed  by  a transformation 
that  carries  points  in  the  local  radar  plane  into  points  on  the  master  plane. 

In  practice,  exact  duplication  of  ideal  conversion  and  transformation  proc- 
esses is  not  possible.  As  a result,  there  are  conversion  errors  and  trans- 
format ior.  errors,  and  these  adversely  affect  the  master  plane  representation 
of  targets. 

Part  I of  the  paper  treats  conversion  error;  i.e.,  the  distance  between  images 
in  the  local  radar  plane  of  target  position  (in  terms  of  altitude,  azimuth, 
and  slant  range)  under  ideal  conversion  and  system  implementations  thereof. 
These  images  are  viewed  as  elements  of  a complex  plane  in  which  the  center  of 
coordinates  represents  the  radar  site.  Thus,  conversion  error  can  be 
expressed  in  terms  of  a range  error  equal  to  the  difference  in  moduli  of  two 
complex  numbers  and  an  angle  error  equivalent  to  the  difference  between  the 
arguments  of  the  same  numbers.  As  the  image  of  target  position  under  ideal 
conversion  traverses  a circle  centered  on  the  origin  of  coordinates,  both 
the  angle  error  and  the  range  error  oscillate  about  median  values.  The  ampli- 

V tude  of  each  oscillation  increases  with  the  radius  of  the  circle.  In  addition, 
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Che  median  range  error  is  strongly  dependent  upon  the  radius  of  the  circle. 

On  the  other  hand,  the  median  angle  error  is  identically  zero. 

In  the  case  where  error  control  is  effected  by  commonly  accepted  procedures, 
the  angle  error  is  less  than  0.0056*  under  standard  operational  conditions. 
Moreover,  the  amplitude  of  the  range  error  is  upper  bounded  by  0.02  nautical 
mile,  whereas  the  median  range  error  can  exceed  0.2  nautical  mile.  However, 
at  least  from  the  theoretical  point  of  view,  there  is  another  method  of 
error  control  capable  of  substantially  eliminating  the  median  range  error 
without  disturbing  the  amplitudes  of  the  angle  and  range  error  oscillations. 
Needless  to  say,  the  reduction  in  error  attainable  with  the  alternative  method 
is  significant.  While  the  practicality  of  the  approach  in  the  present  environ- 
ment is  a moot  question,  the  method  does  offer  some  interesting  possibilities 
in  the  context  of  a future  replacement  for  the  NAS  computer  system. 

Part  II  of  the  paper  deals  with  transformation  error;  i.e.,  the  distance 
between  images  in  the  master  plane  of  the  same  element  in  the  local  radar 
plane  under  ideal  transformation  and  polynomial  approximations  thereof. 

Tight  upper  and  lower  bounds  are  derived  for  the  error  generated  by  the 
nth-order  approximation  to  the  ideal  transformation  equation  under  parameter 
constraints  consistent  with  the  size  of  coverage  areas  associated  with 
current  air  traffic  control  centers  (ARTCC's)  and  the  range  of  operational 
radars.  In  the  case  of  the  first-order  approximation,  the  error  can  exceed 
0.8  nautical  mile,  whereas  the  error  involved  in  the  second-order  approxi- 
mation is  less  than  8 meters.  Clearly,  from  the  standpoint  of  error  control, 
the  first  order  approximation  is  far  less  desirable  than  the  f.econd-order 
approximation.  However,  the  error  introduced  by  the  nth-order  approximation 
is  dependent  upon  the  location  of  the  radar  site  from  which  positional  infor- 
mation is  collected  as  well  as  the  distance  between  the  target  and  the  radar. 
Thus,  as  in  the  current  ATC  environment,  by  using  both  approximations,  it  is 
possible  to  take  advantage  of  the  simplicity  of  the  first-order  approximation 
without  accepting  an  intolerable  transformation  error. 

Although  not  altogether  obvious,  the  transformation  process  is  intimately 
connected  with  the  so-called  magnification  factor.  This  is  just  the  amount 
by  which  the  length  of  an  infinitesimal  arc  on  the  reference  ellipsoid  is 
amplified  under  stereographic  projection  into  the  master  plane.  The  factor 
varies  with  the  position  of  the  arc  on  the  reference  ellipsoid.  For  purposes 
of  convenience,  as  well  as  prediction  accuracy,  it  is  desirable  that  the 
magnification  factor  be  close  to  unity  over  the  coverage  region  of  the  surveil- 
lance system.  In  Part  II,  practical  procedures  for  accomplishing  this 
desideratum  (other  than  those  in  current  use)  are  discussed.  These  procedures 
provide  an  upper  bound  on  the  deviation  of  the  magnification  factor  from 
unity  over  the  coverage  region.  The  bound  is  dependent  on  coverage  region 
size.  For  example,  if  the  maximum  distance  between  any  two  points  in  a coverage 
region  is  on  the  older  of  1,100  nautical  miles,  then  it  is  possible  to  keep 
the  deviation  below  0.004.  As  a result,  the  magnification  factor  must  fall 
between  0.996  and  1.004  over  the  entire  coverage  region.  In  the  case  of 
coverage  regions  with  diameters  less  than  1,100  nautical  miles,  smaller 
bounds  on  the  deviation  from  unity  are  attainable. 
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In  conclusion,  it  is  emphasized  that  this  investigation  has  been  conducted 
under  the  assumption  that  sensor  measurements  are  exact  and  data  processing  is 
carried  out  with  infinite  precision.  In  practice,  this  is  not  the  case. 
Receiver  front-end  noise  exists,  range  and  azimuth  measurements  are  quantized 
before  being  passed  on  to  the  central  processor,  trigonometric  functions  and 
other  standard  mathematical  relationships  are  approximated  in  various  ways  by 
system  hardware  and  software,  and  multiplication  and  addition  are  carried  out 
using  finite  precision  arithmetic.  All  of  these  contribute  to  error,  and  we 
have  no  reason  to  believe  the  effect  upon  conversion  and  transformation  is 
inconsequential.  In  this  sense,  the  results  of  the  paper  should  be  viewed  as 
representative  of  these  processes  under  ideal  operating  conditions. 


1. 


INTRODUCTION 


This  is  the  first  part  of  a two-part  paper  dealing  with  two-dimensional 
representations  of  target  position  from  measurements  of  altitude,  slant  range, 
and  azimuth  in  a multiple  radar  data  processing  system.  It  is  concerned  with 
the  case  in  which  such  representations  are  realized  by  means  of  stereographic 
projection.  Assuming  that  the  mean  sea  level  surface  of  the  earth  is  an 
ellipsoid  of  revolution,  this  technique  can  be  used  to  establish  a unique 
relationship  between  points  on  a plane  and  orthogonal  projections  of  points  in 
space  onto  the  earth's  surface.  The  latter  can  be  determined  from  measurement i 
of  slant  range,  azimuth,  and  altitude  above  mean  sea  level.  Thus,  in  the  case 
of  targets  at  the  same  latitude  and  longitude  but  different  altitudes,  such 
measurements  are  mapped  into  a single  point  on  the  plane.  In  this  way,  data 
obtained  from  many  radars  can  be  used  to  assemble  a single  planar  representa- 
tion of  positions  of  aircraft  and  other  targets  in  three-dimensional  space. 

The  mapping  of  target  positions  onto  a single  plane  can  be  viewed  as  a two- 
stage  procedure.  First,  there  is  a conversion  of  target  altitude  together 
with  target  slant  range  and  azimuth  relative  to  any  given  radar  site  into 
a point  on  a local  radar  plane.  This  is  followed  by  a transformation  that 
carries  points  in  the  local  plane  into  points  on  a single  master  plane. 

Thus,  there  are  as  many  local  planes  as  there  are  radars,  and  each  of  these 
is  mapped  into  the  master  plane  to  establish  the  final  planar  representation 
of  positions  of  targets  within  the  coverage  region  of  the  overall  surveillance 
system. 

Ideally,  conversion  involves  a determination  of  the  orthogonal  projection 
of  target  position  on  the  ellipsoidal  representation  of  the  mean  sea  level 
surface  of  the  earth,  a mapping  of  the  projection  onto  a so-called  conformal 
sphere,  and  stereographic  projection  of  the  spherical  surface  onto  the  local 
radar  plane.  Transformation,  in  the  ideal  sense,  is  realized  by  a bilinear 
transformation  that  takes  the  local  plane  into  the  master  plane.  When  carried 
out  in  this  fashion,  the  process  of  conversion  and  transformation  is  equiv- 
alent to  the  three  steps  of  conversion  alone  with  the  exception  that  the  local 
and  master  planes  are  the  same.  Thus,  under  ideal  conversion  and  transfor- 
mation, the  planar  representation  of  a target  is  invariant  to  the  radar  with 
respect  to  which  slant  range  and  azimuth  are  measured. 

Unfortunately,  in  practice,  it  is  difficult  to  realize  the  ideal  forms  of 
conversion  and  transformation.  As  a result,  modifications  of  these  procedures 
are  encountered  in  operational  systems.  Here  and  in  the  sequel  (refer- 
ence 7),  we  will  investigate  differences  between  images  of  a target  under 
ideal  and  approximate  conversion  procedures  as  well  as  differences  between 
images  of  a point  on  the  local  radar  plane  under  the  ideal  bilinear  transfor- 
mation and  polynomial  approximations  thereof.  While  the  subject  is  far  from 
new  (references  1,  2,  and  3),  numerical  results  are  scarce,  conversion  dif- 
ferences are  expressed  by  formulae  derived  from  spherical  earth  models,  and 
optimization  of  differences  is  defined  in  terms  of  minimax  criteria  applied  to 
targets  above  the  radar  platform  plane.  Objectives  of  this  investigation  are 


1 


to  assemble  a comprehensive  and  yet  concise  picture  of  pertinent  theoretical 
results  applicable  to  an  ellipsoidal  earth  model,  to  provide  numerical  results 
characterizing  differences  between  ideal  procedures  and  those  encountered  in 
practice,  and  to  consider  alternative  methods  for  maintaining  such  differences 
within  tolerable  bounds. 

Our  study  relies  heavily  upon  the  work  of  Flax  (reference  4),  which  is  an 
elegant  unified  treatment  of  theoretical  aspects  of  conversion  and  transforma- 
tion based  upon  the  assumption  of  a spherical  earth.  Here,  his  results  are 
combined  with  the  concept  of  a conformal  sphere  (reference  5)  to  provide 
mathematical  frameworks  for  both  ideal  and  approximate  forms  of  conversion  and 
transformation  in  the  context  of  an  ellipsoidal  earth  model.  In  addition, 
numerical  results  exhibiting  the  differences  between  the  exact  and  approximate 
conversion  procedures  are  examined,  and  attention  is  directed  to  the  problem 
of  reducing  such  differences.  Numerical  investigations  of  transformation  are 
left  to  the  sequel. 

The  next  section  summarizes  essential  features  of  geodetic  coordinates  used 
in  future  developments.  The  relationship  between  geodetic  coordinates  and 
measurements  of  slant  range,  azimuth,  and  altitude  is  discussed  in  section  3. 
Section  4 describes  the  relationship  between  conformal  and  geodetic  latitudes 
and  exhibits  the  mathematical  formulas  by  which  the  surface  of  a conformal 
sphere  can  be  stereographically  projected  onto  a plane.  In  addition,  the 
connection  is  shown  between  stereographic  projection  and  the  planar  repre- 
sentation of  targets  in  multiple  radar  surveillance  systems.  The  details  of 
conversion  are  presented  in  section  5.  In  section  6,  the  differences  between 
the  results  of  ideal  conversion  nnd  those  achieved  by  approximations  employed 
in  operational  systems  are  defined  in  precise  terms.  Section  7 provides 
numerical  results  that  characterize  the  behavior  of  these  differences  for  some 
realistic  target  positions  and  radar  site  locations  in  the  case  where  error 
control  is  effected  by  means  of  a minimax  procedure.  An  alternative  procedure 
for  controlling  errors  is  discussed  in  section  8.  Concluding  remarks  appear  in 
section  9. 


2. 


GEODETIC  COORDINATES 


The  gravity  field  of  the  earth  is  the  result  of  two  forces,  the  centrifugal 
force  due  to  the  earth's  rotation  and  the  force  of  mass  attraction  between  the 
earth  and  other  bodies.  It  can  be  described  in  terms  of  equipotential  sur- 
faces, and  the  force  of  gravity  is  perpendicular  to  such  surfaces.  The 
equipotential  surface  at  mean  sea  level  is  called  the  geoid.  It  is  irregular 
and  hence  unsuitable  as  a basis  for  a coordinate  system.  The  geoid  can  be 
closely  approximated  by  an  ellipsoid  of  revolution  called  the  reference 
ellipsoid.  This  geometry  is  often  used  for  purposes  of  navigation,  gravity 
computations,  etc. 

An  ellipsoid  is  illustrated  in  figure  1.  It  has  a semimajor  axis  of  length  a 
(the  equatorial  radius)  and  a semiminor  axis  of  length  b (the  polar  radius). 

In  the  case  of  the  reference  ellipsoid  (reference  6), 
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FIGURE  1.  REFERENCE  ELLIPSOID 


a = 6,378,388  meters  and  b = 6,356,912  meters. 

The  direction  of  the  earth's  rotation  is  indicated  by  the  arrowhead  (pointing 
north)  on  the  minor  axis  in  conjuntion  with  the  usual  right-hand  rule. 

The  ellipsoid  is  split  into  halves  by  a plane  (the  equatorial  plane)  perpen- 
dicular to  the  minor  axis.  These  correspond  to  the  usual  notions  of  Northern 
and  Southern  Hemispheres.  Planes  parallel  to  the  equatorial  plane  cut  the 
ellipsoid  in  circles  called  parallels.  Planes  containing  the  minor  or  polar 
axis  cut  it  in  ellipses.  The  portion  of  any  such  ellipse  on  either  side  of 
the  polar  axis  is  called  a meridian. 


A point  P in  space  can  be  characterized  by  so-called  geodetic  coordinates  in 
terms  of  meridians,  parallels,  and  the  distance  to  its  perpendicular  projection 
Q on  the  ellipsoid.  Except  in  the  case  of  the  North  or  South  Pole,  the 
projection  is  the  intersection  of  only  one  meridian  and  one  parallel.  In 
the  case  where  the  projection  is  a pole,  the  corresponding  parallel  is  a point 
and  so  there  is  no  ambiguity.  Finally,  the  distance  of  P from  its  projection 
determines  its  location  relative  to  the  parallel-meridian  intersection. 

This  distance,  i.e.,  the  altitude  of  P,  is  regarded  as  negative  when  the  point 
is  enclosed  by  the  ellipsoid  and  positive  otherwise. 


Meridians  and  parallels  are  usually  expressed  as  angles,  namely,  latitude 
and  longitude.  These  can  be  readily  defined  in  terms  of  a Cartesian  system 
with  origin  at  the  center  of  the  ellipsoid.  Referring  to  figure  1,  the 
Cartesian  system  consists  of  the  polar  axis  and  two  orthogonal  axes  in  the 
equitorijil  plane.  ^Positive  directions  along  the  axes  are  indicated  by  unit 
vectors  i,  j,  and  k.  The  meridian  passing  through  the  projection  Q is 
described  by  the  angle  X (longitude)  measured  in  the  equatorial  plane  counter 
to  the  direction  of  the  earth's  rotation  between  the  plane  containing  the 
arrows  designated  as  j and  k and  that  containing  Q and  the  polar  axis. 
Latitude  defines  the  parallel,  and  it  is  the  angle  L in  the  latter  plane 
measured  between  the  equatorial  plane  and  the  line  normal  to  the  reference 
ellipsoid  at  Q.  When  the  projection  lies  in  the  Northern  Hemisphere,  L is 
nonnegative  and  does  not  exceed  90s.  Otherwise,  it  is  negative  but  not  less 
than  -90s. 


As  indicated  in  figure  1,  the  position  vector  of  the  point  P is  just  the  sum 
of  the  position  vector  s(L,X  ) of  its  projection  and  the  outward  unit 
normal  ^(L, X ) to  the  ellipsoid  at  Q multiplied  by  its  altitude  H.  It  can 
be  shown  that 


A A A 

®(L  X )»  a cos  L (sinXi+cosXj)+a(l-e2)sin  L k 
j/l-e2sin2L  / l-e2sin2L 


(1) 


where  e is  the  ellipticity  of  the  ellipsoid;  i.e., 


e2  * 1 - (b/a)2. 


(2) 
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Moreover,  from  figure  1,  it  is  not  difficult  to  see  that 

A A * A 

u3(L,A  ) ■ cos  L (sin  Ai  ♦ cosA  j)  + sin  L k. 


(3) 


Thus,  the  position  vector  of  point  P in  terms  of  its  latitude,  longitude,  and 
altitude  can  be  expressed  as 


s(L, A ) + H U3(L,A  ) 


(4) 


In  later  sections,  we  shall  find  it  convenient  to  refer  to  unit  vectors 
Ql(L,  A)  and  &2(L, A )•  In  figure  1,  these  are  pictured  as  lying  In  the  plane 
tangent  to  the  surface  of  the  reference  ellipsoid  at  Q.  In  addition,  the 
vector  U2(L,A  ) lies  in  the  plane  defined  by  the  A-meridian  and  is  directed 
from  the  South  toward  the  North  Pole.  Inspection  of  figure  1 shows  that 


A A ^ 

U2(L, A ) ■ -sin  L sin  A i -sin  L cos  A j + cos  L k. 
The  remaining  vector  is  defined  by  the  cross  product 
uj(L, A ) * U2(L,  A ) X U3<L, A ) 

* A 

* - cos  A i + sin  X j* 


(5) 


(6) 


It  can  be  viewed  as  lying  in  the  plane  defined  by  the  L-parallel  and  pointing 
along  the  parallel  in  the  direction  of  the  earth's  rotation.  In  the  degen- 
erate case  where  the  projection  of  P is  either  the  North  or  the  South  Pole,  it 
is  sufficient  to  consider  ui(L,A)  and  U2(L,A)  to  be  any  orthogonal 
pair  of  unit  vectors  such  that 


U3(L,A  ) * u} (L, A ) X U2(L,A  ). 


SLANT  RANGE  AND  AZIMUTH 


(7) 


Consider  a radar  at  altitude  Hr  above  the  reference  ellipsoid  at  latitude  Ls 
and  longitude  As.  Let  (L,  A , H)  represent  the  geodetic  coordinates  of  some 
target.  The  position  vector  of  the  radar  site  can  be  expressed  as 

r(Ls,  A.,  Hr)  - s(Ls,  As>  + Hr  «3<L8.  *8),  (8) 

% 

and  the  position  vector  r (L,  A , H)  of  the  target  can  be  written  in  an 
analogous  form.  As  shown  in  figure  2, 

t(L, A,  H)  - r(L,  A,  H)  - r(L8,A8,  Hr)  (9) 

can  be  viewed  as  an  arrow  drawn  from  the  radar  site  to  the  target. 
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Obviously,  the  magnitude  S of  the  vector  t(L,X  , H)  is  just  the  slant 
range  of  the  target.  Ue  can  also  express  target  azimuth  as  a function  of 
this  vector  and  the  unit  vectors  ui(Lg,  X8)  and  u2(L8,  X8).  In  order 
that  the  expression  yield  an  angle  measured  in  the  plane  of  the  radar  platform 
in  accordance  with  the  usual  notions  of  azimuth,  it  is  essential  that  the  axis 
of  rotation  of  the  platform  be  normal  to  the  surface  of  the  ellipsoid  at 
(l>8,  X 8 , 0).  Thus,  in  figure  2,  the  axis  of  rotation  is  depicted  as 
coincident  with  the  unit  normal  vector  U3(L8,  X8).  As  shown  in  the  same 
figure,  it  is  assumed  that  the  platform  rotates  in  the  direction  from 
u2(L8,  X 8 ) to  u}(Ls,  X 8)  in  accordance  with  the  right-hand  rule  used 
in  conjunction  with  the  negative  of  U3(L8,  X8).  With  this  understanding, 
we  define  the  target  azimuth  6 in  terms  of  the  orthogonal  projection  of 
t(L, X , H)  into  the  plane  tangent  to  the  ellipsoid  at  the  orthogonal 
projection  of  the  radar  site,  i.e.,  (L8,  X » 0).  More  precisely,  6 is 
the  angle  measured  in  the  tangent  plane  (wfiich  is  parallel  to  the  radar 
platform  plane)  from  u2(L8,  X8)  in  the  direction  of  platform  rotation  to 
the  projection  of  t(L, X , H).  In  terms  of  vector  dot  products, 

sin0-  x,  H)  • X8)l  (10) 

and 

cose  - i[t(L,X  , H)  • u2(L8,  X8)]  (11) 

where 

/lt(L,  X ,H)  • ^(Ls,  X8)l2+[t(L,  X , H)  • u2(L8,  X 8)]2  <12> 


is  the  magnitude  of  the  projection  of  t(L,  X,  H). 

Theoretically,  it  is  possible  to  determine  target  latitude  and  longitude  from 
slant  range,  azimuth,  and  altitude.  For  example,  one  might  search  through 
latitude-longitude  pairs  (*,  * ) such  that 

| t( - , * ,H)  | - S (13) 

and  pick  one  satisfying  (10)  and  (11).  In  many  situations,  this  is  not  prac- 
tical, and  other  simpler,  but  approximate,  procedures  are  used.  More  will  be 
said  about  the  details  later.  The  point  of  interest  here  is  that  the  objec- 
tive is  to  convert  the  triplet  (S,  6 , H)  into  latitude  and  longitude. 

For  purposes  of  display  and  computational  convenience,  representation  of 
target  position  in  terms  of  latitude  and  longitude  as  a point  on  the  ellip- 
soid is  not  appealing.  Representation  of  position  as  a point  on  some  flat 
surface  is  better.  The  latter  can  be  accomplished  via  a one-to-one  mapping  of 
the  reference  ellipsoid  onto  a plane.  One  such  mapping  can  be  obtained  from 
the  method  of  stereographic  projection  in  conjunction  with  the  concept  of 
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conformal  spheres.  The  baaic  idea  ia  to  map  the  ellipsoid  onto  a sphere  in 
a one-to-one  manner,  and  then  map  the  sphere  onto  a plane  in  the  same  manner. 
The  composition  of  the  two  mappings  is  then  a one-to-one  mapping  of  the 
ellipsoid  onto  the  plane.  The  exact  procedure  is  outlined  in  the  next  section. 


4. 


CONFORMAL  SPHERES  AND  STEREOGRAPHIC  PROJECTION 


Consider  the  ellipsoid  of  figure  1.  For  our  purposes,  any  sphere  with  center 
colocated  with  the  center  of  the  ellipsoid  is  a conformal  sphere.  Since 
the  ellipsoid  is  itself  a sphere  when  the  major  and  minor  axes  have  the 
same  length,  we  are  free  to  talk  about  points  on  a conformal  sphere  in  terms 
of  latitude  and  longitude  in  the  same  way  we  use  these  angles  to  describe 
points  on  the  ellipsoid.  With  this  understanding,  we  are  in  a position  to 
discuss  mappings  that  send  points  on  the  ellipsoid  into  points  on  a sphere. 
One  such  mapping  sends  the  point  Q with  geodetic  coordinates  (L,  A ) on  the 
ellipsoid  into  the  point  (<f>,  \ ) on  a conformal  sphere  via  the  relationship 
(references  2 and  5) 


The  supping  is  one  to  one.  Moreover,  Q and  its  image  have  the  same  longitude 
X but  different  latitudes.  The  latitude  <f>  of  the  image  is  often  referred 
to  as  the  conformal  latitude  of  Q to  distinguish  it  from  the  geodetic 
latitude  L. 

As  already  pointed  out,  a sphere  is  a special  case  of  an  ellipsoid.  Hence, 
the  notation  developea  earlier  in  conjunction  with  figure  1 is  applicable  to 
discussions  of  conformal  spheres.  Accordingly,  U3(<^,  A 8 ) can  be  viewed 
as  a unit  vector  normal  to  the  surface  of  the^sphere  illustrated  in  figure  3 
at  latitude  0a  and  longitude  X8.  Likewise,  u^(08,  As)  and  U2 (<£s,  X8) 
can  be  represented  as  arrows  lying  in  a plane  perpendicular  to  the  surface 
normal  at  an  altitude  hr  above  the  spherical  shell.  Adopting  the  convention 
that  U2 X8)  indicates  the  direction  of  increasing  imaginary  parts  of 
complex  numbers,  chis  representation  defines  a complex  plane.  In  the  follow- 
ing paragraphs,  we  will  outline  the  method  of  stereographic  projection  as  a 
mapping  from  the  spherical  surface  onto  this  plane. 

The  line  normal  to  the  sphere  at  X8,  0)  passes  through  the  center 

and  intersects  the  spherical  surface  again  at  (-^fe,  X8  + 1t,  0).  The  latter 
is  the  focal  point  for  stereographic  projection.  Under  the  projection  map- 
ping, the  image  of  a point  A on  the  spherical  shell  is  just  the  point  I where 
the  plane  intersects  the  line  through  A and  the  focal  point. 
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FIGURE  3 


STEREOGRAPHIC  PROJECTION 


Suppose  <f>  is  Che  latitude  and  X the  longitude  of  A,  and  x and  y represent 
the  real  and  imaginary  parts  of  the  complex  representation  of  I.  Also,  for  the 
ament,  suppose  that  the  projection  plane  is  tangent  to  the  sphere;  i.e.,  hr 
is  0.  Then,  assuming  the  radius  of  the  sphere  is  Es,  it  can  be  shown  that 
stereographic  projection  is  equivalent  to  the  mapping  (reference  4) 

(0,  X )+*+iy  (15) 

where 

2ES  cos  <f>  8 in  (X8-X)  (16) 

X l+s in  (f>B  sin  0 +cos  <f>s  cos  <p  cos  (X8-X)  » 

2ES  [ cos  <pa  sin0-sin  <ps  cos  0cos(X8-X) ] (17) 

l+s  in  <PB  sin  <j>  + cos  0s  cos  <p  cos(X8-X) 

This  is  a one-to-one  relation.  The  inverse  mapping  is 


and  i is  ">J-i  . 

defined  by 


sin0  ■ 


(4E2-x2-y2)sin$+4Esy  cos 
4E2+x2+y2 


(18) 


sin  (X-XB )■ 


-4E»x 


(19) 


\/ (-4E8x]2+[ (4E2-x2-y2)  cosqij-XEg  y sim^]2 


and 


cos  (X-Xg)  “ 


(4E2-x2-y2)  cos08-4E8y  sin  <pe 


\A(-4E8x}2+[(4eJ  -x2-y2)  cos08-4EBy  sin08]2 


(20) 


The  preceding  equations  characterise  stereographic  projection  in  the  case 
where  the  projection  plane  is  tangent  to  the  sphere.  Assume  now  that  it  is 
elevated  at  altitude  hr.  Then  the  image  of  A is  just  a scalar  multiple  of 
the  complex  number  that  is  obtained  when  the  altitude  vanishes.  The  multi- 
plicative factor  is 

1+^r  . (21) 

2E, 


In  the  case  of  the  inverse  procedure,  this  factor  is  divided  into  the  complex 
number  obtained  from  the  elevated  plane,  and  then  the  equations  (18)-(20)  are 
applied  to  the  result. 

Ue  now  return  to  the  case  where  the  projection  plane  altitude  hr  is  0.  Let 
us  refer  to  the  complex  plane  tangent  to  the  conformal  sphere  of  radius  Es  at 
latitude  <f>a  and  longitude  Aa  as  the  z-plane.  In  addition,  let  us  consider 
another  conformal  sphere  of  radius  Er  and  a complex  plane  tangent  to  it 
at  latitude  <f>r  and  longitude  Ar.  Ue  will  call  the  latter  plane  the 
w-plane.  It  has  been  shown  that  the  mapping  defined  by  (15)-(17)  stereo- 
graphically  projects  longitude  and  conformal  latitude  into  the  z-plane. 
Likewise,  by  making  the  obvious  substitutions,  the  same  equations  can  be  used 
to  map  longitude  and  conformal  latitude  into  the  w-plane.  Since  both  mappings 
are  one-to-one,  an  obvious  composition  of  the  second  and  the  inverse  of  the 
first  defines  a one-to-one  mapping  of  the  z-plane  onto  the  w-plane.  Under 
this  mapping  and  subject  to  the  condition  that 

|\r-\8|<90\  (22) 

it  can  be  shown  that  the  image  w of  an  element  z is  given  by  the  formula 
(reference  A) 


_JL  z e +wQ 
w " £s 

1 - z e “*®w£ 

*SrBs 

where  w*  is  the  conjugate  of  wQ,  |B|  does  not  exceed  90°,  and 
tan  8 * - sin  (Ar-As)  (sin</y*-sin<£8> 


(23) 


( l+s  in  <f>  rs  in^  )cos  (Ar-As  )+coB<f*ecoa<f>a 


(24) 


Obviously,  Che  complex  number  w0  is  just  Che  image  of  Che  origin  of  che 
z-plane.  However,  ic  is  clear  from  Che  definicion  of  Che  mapping  ChaC  any 
elemenc  of  Che  z-plane  and  ics  image  represenC  Che  same  longiCude  and  laCiCude. 
ConsequenCly,  Che  real  and  imaginary  pares  of  vQ  can  be  found  from  (16)  and 
(17)  by  replacing  Che  subscripC  a wich  r.  As  will  be  seen  shorcly,  Che  con- 
sCrainC  (22)  is  noC  difficulC  Co  meeC  in  many  pracCical  siCuaCions. 

EquaCion  (23)  is  imporCanC  in  some  applicaCions  where  measuremenCs  Caken  from 
a mulcicude  of  far-flung  radars  are  used  for  surveillance  purposes  aC  a single 
site.  As  will  be  seen  in  Che  following  secCion,  estimates  of  longiCude  and 
conformal  laCiCude  from  measuremenCs  of  slanC  range  and  azimuch  by  a single 
radar  in  conjuncCion  wich  reporCs  of  CargeC  alCiCude  can  be  formed  by  means  of 
a mapping  inCo  a complex  plane  CangenC  Co  an  appropriaCe  sphere  ac  Che  longi- 
Cude and  conformal  laCiCude  of  Che  radar  sice.  This  plane  corresponds  Co  Che 
z-plane,  and  Chere  are  as  many  such  complex  planes  as  Chere  are  radars.  Each 
of  Chese  is  mapped  inCo  a single  masCer  plane,  corresponding  Co  Che  w-plane, 
via  a funcCional  relationship  analogous  Co  (23).  Thus,  images  on  Che  masCer 
plane  represenC  longicudes  and  conformal  laCiCudes  of  CargeCs  wichin  Che  com- 
pos ice  of  Che  coverage  regions  of  many  individual  radars.  Clearly,  Che  con- 
aCrainC  (22)  will  be  sacisfied  so  long  as  Che  absoluCe  difference  beCween  Che 
longicudes  of  poinCs  of  cangency  of  Che  masCer  plane  and  Che  local  radar  plane 
does  noC  exceed  90*  for  every  radar. 


CONVERSION 


In  Chia  secCion,  we  consider  a radar  aC  longiCude  As,  geodeCic  latitude 
Ls,  and  alCiCude  Hr  relative  Co  Che  reference  ellipsoid.  We  also  consider 
a target  at  altitude  H for  which  Che  slant  range  and  azimuth  relative  to  Che 
radar  are  S and  e , respectively.  As  pointed  out  earlier,  CargeC  latitude 
and  longiCude  can  be  determined  from  S,  6,  and  H.  Here,  we  review  how  this 
is  accomplished  in  some  practical  situations  by  means  of  an  approximation 
based  upon  Che  method  of  stereographic  projection  and  Che  assumption  of  a 
spherical  earth. 

In  the  vicinity  of  the  radar,  the  surface  of  the  reference  ellipsoid  is  some- 
what like  ChaC  of  a sphere  of  radius 


E8-  | s(L8,A8)|  , (25) 

the  distance  between  the  center  of  the  ellipsoid  and  the  orthogonal  projection 
of  Che  radar  site  onto  Che  ellipsoidal  surface.  The  geometry  is  illustrated  in 
figure  4 where  Che  radar  is  pictured  as  a point  removed  a radial  distance  Hr 
from  the  spherical  shell.  This  configuration  is  used  Co  develop  equations 
employed  in  radar  systems  for  converting  slanC  range,  azimuth,  and  altitude 
into  latitude  and  longitude  or  the  equivalent  thereof. 

Consistent  with  the  spherical  geometry  of  figure  4,  the  orthogonal  projection 
of  the  target  T onto  the  geoid  is  represented  by  the  intersection  A of  the 
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spherical  shell  and  the  radial  line  passing  through  the  target.  Points  on  the 
shell  are  stereographically  projected  into  a plane  perpendicular  to  the  axis 
of  rotation  of  the  radar  platform  at  altitude  hr.  This  is  accomplished  by 
lines  through  the  intersection  F of  the  axis  and  the  shell  furthest  removed 
from  the  radar  site.  The  image  of  A under  the  projection  mapping  is  repre- 
sented by  the  point  I,  the  intersection  of  the  plane  and  the  line  through  F 
and  A. 

The  location  of  I in  the  projection  plane  can  be  expressed  in  terms  of  the 
coordinate  system  described  in  the  preceding  section.  Letting 


R " Vs2-(h-Hr)2 


(26) 


it  can  be  shown  that  the  coordinates  of  I are  just  the  real  and  imaginary 
parts  of  the  complex  image  of  slant  range,  azimuth,  and  altitude  under  the 
mapping 

(S,  d , H)  + De  (27) 

where  (reference  4) 


(i+hr  ) 

2ES  K • (28) 


\/(i+|  ) dA)V 

8 Es  4E2 
s 


Obviously,  the  coordinates  of  I can  be  used  to  find  the  latitude  and  longitude 
of  the  target  and  vice  versa  via  equations  (15)-(21).  However,  the  spherical 
viewpoint  ignores  the  difference  between  geodetic  and  conformal  latitudes,  and 
so  we  must  expect  some  error  within  the  context  of  a theory  based  upon  an 
ellipsoidal  representation  of  the  geoid.  Nevertheless,  we  do  expect  it  to 
provide  reasonable  estimates  of  target  longitude  and  conformal  latitude  in  the 
case  where  the  target  is  close  to  the  axis  of  rotation  of  the  radar. 


Actual  implementation  of  the  mapping  (27)  may  not  be  feasible  due  to  limita- 
tions imposed  by  other  than  theoretical  considerations.  In  such  cases,  the 
mapping 

(S,  e , H)+  RE  (29) 

is  often  a satisfactory  substitute.  The  accuracy  of  this  formula  relative  to 
what  can  be  achieved  with  the  mapping  (27)  depends  upon  the  altitude  hr  of 
the  projection  plane.  This  is  readily  apparent  from  figure  4 which  exhibits 
the  separation  of  the  complex  numbers 


Dei(r0) 


and  Rei(r*> 


(30) 


14 


in  Che  plane  as  well  as  Che  reladve  posicion  of  Che  corresponding  images  A 
and  B on  Che  spherical  surface  under  Che  projeccion  mapping.  Happily,  if  we 
accepC  Che  exisCence  of  some  maximum  alcicude  HM  for  CargeCs  of  inCeresC, 

Chen  ic  is  possible  Co  find  a projeccion  plane  alcicude  ChaC  is  opcimal  in  a 
minimax  sense. 

In  parCicular,  referring  Co  figure  4,  ic  can  be  shown  ChaC  R increases  mono- 
Conically  from  a minimum  value  Rq,  to  a maximum  value  Rjj  as  Che  CargeC  T 
moves  ouCward  along  a radial  line  from  Che  radar  plaCform  plane  Co  Che  maximum 
alCiCude  H^.  These  exCreme  values  of  R are  funcCions  of  Che  angle  a beCween 
Che  radial  line  Chrough  Che  CargeC  and  Che  axis  of  rocacion  of  Che  radar. 

When 


2 l (31) 

H 3 H H H H "3 

cos  a m ( 1 ♦ _R)  ( 1 ♦ M ♦ R + *LJL_) 

E_  E.  „2 


Che  difference 

Rm-R®  <32> 

is  maximal  (reference  4).  From  figure  4,  we  expecC  ChaC  chere  exists  a pro- 
jeccion plane  alcicude  for  which  Che  image  I of  a CargeC  on  Che  maximizing 
radial  spliCs  Che  difference  into  equal  pares.  Indeed,  this  is  Che  case,  and 
Che  corresponding  alcicude  is  (reference  4) 

1 23  (33) 

hr«2Es {| J(1+hR+hM+hRhM)3  ♦ 1(1*hR)312_1} 


In  face,  for  any  given  projeccion  plane,  | D-R  | assumes  a maximum  over  all  tar 
get  positions  below  altitude  and  above  Che  radar  platform  plane,  and  it 
can  be  shown  that  the  maximum  is  least  when  the  projection  plane  altitude  is 
given  by  (33). 

The  observations  of  the  preceding  paragraphs  can  be  used  to  show  that  the 
altitude  (33)  is  near  optimal  in  another  sense.  Consider  the  angle  C in 
figure  4 between  the  lines  that  intersect  in  the  stereographic  projection 
center  F and  cut  the  complex  plane  at  the  points  represented  by  the  numbers 
(30).  At  a give  projection  plane  altitude,  this  angle  assumes  a maximum 
value  over  Che  CargeC  positions  below  altitude  HM  and  above  Che  radar  plat- 
from  plane.  On  the  other  hand,  the  angle  of  the  radial  for  which  the  dif- 
ference between  R,,,  and  R^  is  maximal  is  small  when  the  right  side  of  (31) 
is  close  to  1.  In  cases  such  as  this,  it  is  not  difficult  to  argUe  that  the 
maximum  value  of  £ is  essentially  minimized  by  the  projection  plane  at  the 
altitude  (33). 
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6. 


CONVERSION  ERRORS  - DEFINITIONS 


Consider  a target  at  altitude  H,  longitude  X , and  geodetic  latitude  L relative 
to  the  reference  eilipaoid.  Let  S and  6 be  the  slant  range  and  azimuth, 
respectively,  of  the  target  relative  to  the  radar  described  in  the  previous 
section  at  altitude  Hr,  longitude  \s,  and  geodetic  latitude  Ls.  Also,  let 
us  agree  to  use  "z-plane"  to  denote  the  complex  plane  tangent  to  the  conformal 
sphere  of  radius  Es,  defined  by  (25),  at  the  longitude  and  conformal  lati- 
tude <f> 8 of  the  radar.  Finally,  suppose  the  conversion  process  is  accom- 
plished by  mapping  (S,  6,  H)  into  the  point 


R£  (34) 

2E8 


of  the  z-plane  where  R is  defined  by  (26). 

Let  represent  the  conformal  latitude  of  the  target  and 

z«x+iy  (35) 

be  the  z-plane  image  of  the  point  (.<p,  A ) on  the  conformal  sphere  under  the 
mapping  (15).  From  (14)-(17),  it  is  clear  that  the  relation 

(L,  A )-*z  (36) 

defines  a one-to-one  mapping  M from  the  ellipsoid  onto  the  z-plane.  The  map- 
ping is  essentially  a transformation  of  geodetic  latitude  into  conformal 
latitude  followed  by  stereographic  projection  of  the  conformal  sphere  onto  the 
z-plane.  Obviously, 

M(LS,  A 8)“0,  (37) 

since  the  longitude  and  geodetic  latitude  of  the  radar  site  must  map  into  the 
origin  of  the  z-plane. 
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Since  M is  one-to-one,  the  point  z uniquely  determines  target  longitude  and 
geodetic  latitude.  Moreover,  these  angles  together  with  target  altitude 
can  be  used  to  find  target  slant  range  and  azimuth  via  (8)-(12).  Finally,  the 
image  zc  of  the  target  in  the  z-plane  under  the  conversion  process  (34)  is 
determined  by  target  altitude,  slant  range,  and  azimuth.  Thus,  for  a given 
radar  configuration  defined  in  terms  of  maximum  target  altitude  and  the 
geodetic  coordinates  of  the  site  itself,  zc  can  be  viewed  as  a function  of 
target  altitude  and  the  image  z of  the  target  in  the  z-plane  under  the  mapping 
M.  Indeed,  this  is  the  viewpoint  that  is  adopted  in  the  following  paragraphs. 

Ideally,  the  objective  of  conversion  is  to  map  (S,  6,  H)  into  z.  However,  due 
to  geometric  and  analytical  approximations,  the  complex  numbers  z and  zc  are 
not  necessarily  the  same.  The  difference 

6 j(z,H)=  |z-zc  | , (38) 

for  want  of  a better  name,  will  be  called  the  conversion  error.  As  indicated 
in  figure  5,  it  is  also  convenient  to  talk  about  the  errors 

62(z,H)=|z  |-  |zc|  (39) 

and 

53(z ,H)=arg(z)-arg(zc)  (40) 

where  "arg"  denotes  the  argument  of  a complex  number. 


We  will  have  occasion  to  use  the  polar  forms  of  the  z-plane  images  of  the 
target  under  the  mapping  M and  the  conversion  process,  namely 


z*  pe 


iy 


and 


(41) 


We  will  call  p and  Y the  rang*  and  angle  of  the  target,  whereas  pc  and  y c 
will  be  referred  to  as  the  rsnge  and  angle  due  to  conversion.  Thus,  (39)  and 
(40)  can  be  viewed  as  errors  in  range  and  angle  that  together  give  rise  to  the 
overall  coversion  error  (38).  In  what  follows,  we  will  often  aMse  the  nota- 
tion introduced  in  (38)-(40)  by  deferring  to  the  polar  form  of  complex  num- 
bers and  writing  ^(p,  Y , H)  in  place  of  ^(z,  H)  for  k=l,  2,  or  3. 


Let  us  now  consider  the  situation  in  which  target  altitude  remains  constant 
and  the  trajectory  induced  by  target  motion  in  the  z-plane  under  the  mapping 
M is  a circle  centered  at  the  origin.  This  is  equivalent  to  holding  p and  H 
constant  while  allowing  Y to  traverse  the  continuum  between  0 and  2 it  radians. 

In  general,  the  error  ‘^(p,  Y»  H)  will  oscillate  about  the  median 

m^PjH)®  [max  6]c(p,Y,H)+min  5 ^(  p,  Y,H)  )/2  (42) 

0<Y<2u  0<y<2tt 
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(43) 


with  an  amplitude* 


ajc(p,H)*max  ,Y,H)-nnc(p,H) 

0<  y <2v 

For  example,  figure  6 illustrates  the  behavior  of  the  range  and  angle  errors 
when  target  altitude  is  40,000  feet,  the  target  range  is  160  nautical  miles 
(i  nautical  mile  ■ 1,832  meters)  and  hr  satisfies  (33).  As  shown,  the 
median  range  error  is  .0047  nautical  mile  accompanied  by  an  amplitude  of  .012 
nautical  mile. 


7. 


CONVERSION  ERRORS  - NUMERICAL  RESULTS 


In  this  section,  we  will  consider  some  numerical  results  obtained  for  conver- 
sion errors  in  connection  with  radar  site  locations  in  the  Northern  Hemisphere 
on  parallels  28,  38,  and  48  at  altitudes  of  0 and  10,000  feet.  In  each  case, 
the  maximum  altitude  for  targets  of  interest  is  60,000  feet.  Also,  the 
elevation  of  the  projection  plane  is  given  by  the  right  side  of  (33),  and  so 
it  tends  to  minimize  the  maximum  conversion  error  over  targets  above  the  radar 
platform  plane  and  below  H^.  However,  it  is  well  to  recognize  that  little 
is  to  be  gained  by  examining  errors  associated  with  targets  that  are  hidden 
from  the  field  of  view  of  the  radar.  Thus,  rather  than  proceeding  directly 
with  a discussion  of  errors,  something  will  be  said  first  about  what  we  mean 
by  radar  horizon. 

Actual  radars  are  capable  of  detecting  targets  somewhat  below  the  line-of- 
sight  horizon.  Nevertheless,  it  is  an  obvious  consequence  of  the  mathematical 
setup  used  in  this  paper  that  a necessary  condition  for  a target  to  be  visible 
to  the  radar  is  that  the  line  segment  connecting  the  two  entities  does  not 
intersect  the  reference  ellipsoid  in  more  than  one  point.  More  specifically, 
consider  the  plane  defined  by  the  axis  of  rotation  of  the  radar  platform 
and  a line  through  the  radar  tangent  to  the  ellipsoid.  The  line  consists  of 
two  semi-inf inite  segments  terminating  in  the  point  of  tangency.  One  of  these 
contains  the  radar  site.  As  shown  in  figure  7,  targets  in  the  shaded  region 
of  the  plane  below  the  other  segment  and  above  the  surface  of  the  ellipsoid 
cannot  be  seen  from  the  radar  site. 

The  collection  of  semi-infinite  line  segments  that  are  tangent  to  the  ellip- 
soid and  emanate  from  the  radar  site  form  a cone-like  surface.  We  will  refer 
to  the  locus  of  points  at  which  it  intersects  the  ellipsoid  as  the  radar 
horizon  contour.  The  contour  divides  the  cone  into  a surface  of  infinite 
extent  and  one  of  finite  extent.  We  will  refer  to  the  union  of  the  surface  of 
infinite  extent  with  the  orthogonal  projection  of  the  surface  of  finite  extent 
onto  the  ellipsoid  as  the  radar  horizon  surface.  As  indicated  in  figure  7,  a 
target  above  the  reference  ellipsoid  can  be  seen  from  the  radar  site  only  if 
it  is  also  above  the  horizon  surface. 


*In  the  case  of  a spherical  earth  model  the  amplitude  a^  (p,  H)  vanishes 
identically  for  all  k. 


19 


RANGE  ERROR  5 ( p 


FIGURE  6.  RANGE  AND  ANGLE  ERRORS  FOR  RADAR  AT  SEA  LEVEL  ON  PARALLEL  38 


As  the  point  z moves  away  from  the  origin  of  the  z-plane  along  a radial  path, 
we  expect  the  altitude  of  the  horizon  surface  at  the  longitude  and  geodetic 
latitude  represented  by  M~l(z)  to  be  a non-decreasing  function  of  p.  In 
fact,  we  expect  the  function  to  differ  for  different  radial  paths.  Thus,  in 
order  to  characterize  the  horizon  surface  by  means  of  a single  conservative 
descriptor,  it  is  natural  to  ask  for  a radial  path  along  which  the  horizon 
altitude  increases  at  the  greatest  rate.  For  our  purposes,  we  are  Interested 
in  such  optimal  paths  in  the  case  where  the  radar  horizon  contour  lies  wholly 
in  the  Northern  Hemisphere.  Unfortunately,  at  least  to  within  the  limits  of 
our  knowledge,  this  is  an  unsolved  problem.  On  the  other  hand,  a cursory 
examination  of  simple  examples  leads  us  to  believe  that  the  rate  of  increase 
is  greatest  along  the  negative  portion  of  the  Imaginary  axis.  However,  this 
is  a matter  of  conjecture  rather  than  proof.  Nevertheless,  right  or  wrong, 
the  choice  provides  all  we  need,  namely,  a rough  idea  of  what  targets  fall 
within  the  field  of  view  of  the  radar.  The  horizon  surface  altitude  function 
corresponding  to  the  selected  z-plane  path  is  exhibited  in  figure  8 for  a 
radar  at  sea  level  on  parallel  38.  Figure  9 depicts  the  function  in  the 
case  where  the  radar  site  elevation  is  10,000  feet  on  the  same  parallel. 
Corresponding  curves  representing  horizon  surface  altitudes  for  radars  on 
parallels  28  and  48  are  little  more  than  reproductions  of  figures  9 and  10. 
This,  of  course,  is  to  be  expected  in  view  of  the  near  spherical  geometry  of 
the  reference  ellipsoid. 

Returning  to  the  subject  of  conversion  errors,  consider  two  targets  at  the 
same  altitude  such  that  the  z-plane  representation  of  one  is  the  negative  of 
the  complex  conjugate  of  the  other.  An  obvious  consequence  of  the  conversion 
process  is  that  the  corresponding  angle  errors  are  equal  in  magnitude  and 
opposite  in  sign.  Thus,  the  median  of  any  angle  error  oscillation  vanishes 
identically  as  indicated  by  the  example  of  figure  6. 

Another  consequence  of  the  conversion  process  is  that  angle  error  is  indepen- 
dent of  the  radar  site  altitude  Hr  and  the  projection  plane  height  parameter 
hr.  On  the  other  hand,  it  does  depend  upon  target  altitude  and  the  target 
range.  In  figure  10,  amplitudes  of  angle  error  oscillations  associated  with 
a target  altitude  of  60,000  feet  are  plotted  against  the  target  range  for 
radar  sites  located  on  parallels  28  and  48.  Extensive  examination  of  numeri- 
cal results  for  target  altitudes  between  1,000  and  60,000  feet  indicates  that 
angle  error  amplitude  is  essentially  a nondecreasing  function  of  altitude 
when  the  range  is  fixed.  Consequently,  we  have  treated  figure  10  as  an 
illustration  of  upper  bounds  for  angle  errors. 

Amplitudes  of  range  error  oscillations  are  plotted  against  range  in  figure  11 
for  radar  sites  located  on  parallels  28  and  48.  The  illustration  implies  that 
range  error  amplitude  is  invariant  to  target  altitude  and  the  elevation  of  the 
radar  site.  Strictly  speaking,  this  is  not  so.  However,  numerical  results 
obtained  for  radars  at  0 and  10,000  feet  and  target  altitudes  between  1,000 
and  60,000  feet  indicate  the  variation  is  less  than  a few  tenths  of  one  meter. 
For  our  purposes,  such  minute  differences  are  irrelevant. 


ALTITUDE  IN  FT 


FIGURE  9. 


ALTITUDE  OF  HORIZON  SURFACE  FOR  RADAR  AT  10,000  FEET  ON  PARALLEL  38 


FIGURE  10.  UPPER  BOUNDS  ON  ABSOLUTE  ANGLE  ERROR 


FIGURE  II.  RANGE  ERROR  AMPLITUDE 


Figures  12  and  13  illustrate  medians  of  range  error  oscillations  at  several 
different  target  altitudes  for  radars  at  sea  level  and  10,000  feet.  The  data 
were  obtained  for  radar  site  locations  on  parallel  38.  However,  results 
obtained  for  parallels  28  and  48  are  the  same  to  within  .001  nautical  mile. 
Consequently,  for  our  purposes,  the  illustrations  are  equally  applicable  to 
all  three  parallels.  Actually,  this  Is  not  surprising.  In  the  case  of  a 
spherical  earth  model,  error  is  invariant  to  radar  site  latitude,  and  we 
expect  this  invariance  to  carry  over  in  some  approximate  sense  to  the  model 
embodied  in  the  reference  ellipsoid. 

Figures  14  and  15  illustrate  the  maximum  overall  conversion  error  for  radars 
on  parallel  36.  As  already  pointed  out,  the  projection  plane  elevation  used 
in  this  section  can  be  regarded  as  a solution  to  a minlmax  problem  involving 
targets  in  the  region  above  the  radar  platform  plane  and  below  60,000  feet. 
This  region  corresponds  to  that  portion  of  figures  14  and  15  below  and  to 
the  right  of  the  curve  associated  with  targets  at  60,000  feet.  Thus,  the 
optimality,  if  any,  of  the  data  above  and  to  the  left  of  this  boundary  is  an 
open  question.  When  figures  14  and  15  are  compared  with  figures  12  and  13, 
it  becomes  apparent  that  median  range  error  contributes  heavily  to  the  over- 
all conversion  error.  In  the  next  section,  we  will  discuss  a method  for 
reducing  the  size  of  the  range  error. 


8. 


ALTERNATIVE  APPROACHES  TO  ERROR  CONTROL 


The  surface  of  the  reference  ellipsoid  in  the  vicinity  of  the  radar  is  closely 
approximated  by  that  of  a sphere  with  radius  Es  defined  by  (25).  If,  in  our 
mind's  eye,  we  allow  the  lengths  of  the  semiminor  and  semimajor  axes  of  the 
ellipsoid  to  approach  E8,  then  angle  errors  and  amplitudes  of  range  error 
oscillations  vanish.  All  that  is  left  is  the  median  range  error  m2<p,  H). 

In  fact  z and  zc  of  (41)  can  be  obtained  by  multiplying  the  complex  numbers 
(30)  by  the  Inverse  of  1 + hr/2Es.  Thus,  in  the  context  of  a spherical 
representation  of  the  geoid  of  radius  Es, 

m2(p,  H)  - p - R [1  + hr/2Es]-l  (44) 


where 


P - D [I  + hr/2E8]"1 


(45) 


and  R is  related  to  D through  the  expression  (28).  From  (28)  and  (45),  R can 
be  expressed  as  a function  of  p , and  this  can  be  used  to  reduce  (44)  to  the 
form 


m2(p,  H) 


V(1  + H/Eb)  (1  + Hr/Eb) 

(1  + hr/2Es)  Vl  + p2/4eJ 


(46) 
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FIGURE  14.  MAXIMUM  CONVERSION  ERROR  FOR  RADAR  AT  SEA  LEVEL  ON  PARALLEL  38 
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FIGURE  15.  MAXIMUM  CONVERSION  ERROR  FOR  RADAR  AT  10,000  FEET  ON  PARALLEL  38 
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As  indicated  in  Section  5,  the  right  side  of  (46)  is  optimized  in  a minimax 
sense  when  hr  satisfies  (33).  In  addition,  it  can  be  used  to  approximate 
the  median  range  error  data  in  figures  12  and  13  derived  from  the  setup  of 
Section  6 based  on  the  reference  ellipsoid.  Moreover,  the  approximation  is 
essentially  invariant  to  Eg.  In  particular,  using  any  value  for  Es  between 
the  polar  radius  and  the  equatorial  radius  of  the  reference  ellipsoid,  dif- 
ferences between  the  median  range  errors  generated  by  (46)  and  the  ellipsoidal 
earth  model  are  hardly  discernable  to  the  naked  eye  when  the  errors  are 
exhibited  as  a graph  with  a scale  identical  to  that  adopted  in  figures  12 
and  13.  From  this  viewpoint,  it  is  apparent  that  median  range  error  is  con- 
trolled in  a minimax  sense  when  hr  is  selected  in  accord  with  the 
formula  (33). 


As  will  be  seen,  there  are  other  ways  to  control  error  for  which  hr  does  not 
satisfy  (33).  To  this  end,  it  is  important  to  recognize  several  characteris- 
tics of  the  relationship  (46).  First,  if  hr  is  determined  from  the 
expression 


1 + hr/2E8 


+ K/E8)  (1  + HR/E8) 


(47) 


where  K is  some  constant,  then  the  right  side  of  (46)  is  invariant  to  the 
altitude  of  the  radar  site.  Second,  if  our  observations  of  the  preceding 
paragraph  regarding  invariance  to  Es  also  apply  to  this  alternative  method 
of  selecting  the  projection  plane  altitude  hr,  then  the  right  side  of  (46) 
is  essentially  invariant  to  all  aspects  of  the  location  of  the  radar  site. 
Third,  if  in  addition,  the  right  side  of  (46)  closely  approximates  the  median 
range  error  associated  with  the  new  selection  process  for  hr  in  the  context 
of  the  setup  of  Section  6,  then  we  will  have  succeeded  in  achieving  a median 
range  error  that  is  essentially  invariant  to  radar  site  location.  In  fact, 
the  median  range  error  will  then  be  amply  represented  by  any  one  of  the 
numbers 


f (p,  H,  E) 


VI  ♦ H/E 

1 + K/E 


Vl  ♦ p2/4E2 


where  E is  not  less  then  the  polar  radius  nor  greater  than  the  equatorial 
radius  of  the  reference  ellipsoid,  i.e., 


(48) 


b < E < a 


(49) 


If  all  of  this  is  fact,  then,  by  replacing  P with  the  apparent  range  pc  observed 
as  a result  of  the  conversion  process,  the  concept  of  adding  f (pc,  H,  E)  to 
pc  for  the  purpose  of  correcting  median  range  error  has  some  rather  obvious 
practical  implications  as  a means  for  controlling  conversion  error. 


Returning  now  to  the  setup  of  section  6,  let  us  consider  the  situation  where 

K - 61,000  feet  (50) 


and  Che  projection  plane  altitude  satisfies  (47).  Figure  16  depicts  the 
median  range  error  as  a function  of  target  range  and  altitude  when  the  radar 
site  is  at  mean  sea  level  on  parallel  38.  We  have  purposely  neglected  to 
exhibit  the  location  of  the  radar  in  the  illustration.  The  reason  is  that 
when  similar  curves  are  plotted  on  the  same  graph  paper  for  other  radar  site 
locations,  differences  are  not  discernible  to  the  naked  eye.  In  particular, 
only  differences  less  than  .001  nautical  mile  were  observed  among  radar  site 
locations  at  all  six  combinations  of  parallels  28,  38,  and  48  and  altitudes  of 
0 and  10,000  feet.  Moreover,  in  all  these  cases,  plots  of  amplitudes  of  target 
range  oscillations  appear  as  copies  of  figure  11,  and  actual  numerical  dif- 
ferences are  less  than  a few  tenths  of  one  meter.  In  addition,  as  pointed  out 
earlier,  angle  error  is  strictly  invariant  to  the  projection  altitude  hr  so 
that  figure  10  represents  the  angle  errors  in  the  immediate  situation.  Thus, 
except  for  some  rather  minute  perturbations,  the  sole  result  of  choosing  the 
projection  plane  altitude  in  accordance  with  (47)  is  the  creation  of  a func- 
tional relationship  between  median  range  error,  target  range,  and  target 
altitude  that  is  invariant  to  the  latitude  and  altitude  of  the  radar  site. 
Moreover,  it  can  be  readily  verified  that  the  functional  relationship  is 
extremely  well  represented  by  the  expression  (48)  for  any  constant  value  of 
E satisfying  the  constraint  (49).  In  other  words,  the  desiderata  of  the  pre- 
ceding paragraph  have  been  realized,  and,  in  addition,  it  is  clear  that 
figures  10  and  11  represent  upper  bounds  on  the  residual  error  after  the 
median  range  error  has  been  removed. 

Turning  to  figure  16,  it  is  apparent  that  | p — Pci  be  a good  deal  less 

than  1 nautical  mile  under  rather  practical  limitations  on  target  range  p and 
target  altitude  H.  From  the  same  figure,  it  appears  that  the  corresponding 
difference  between  f(p,  H,  E)  and  f(pc,  H,  E)  is  quite  small  compared  to 
f(p,  H,  E).  The  obvious  conclusion  is  that  median  range  error  can  be  essen- 
tially eliminated  by  letting 

zc  - [pc  + f(pc,  H,  E)]e  l7c  (51) 

represent  target  position  in  the  z-plane  instead  of  zc  (see  (34)  and  (41)). 
Moreover,  we  expect  that  a2(p,  H)  of  figure  11  is  very  close  to  being  an 
upper  bound  on  the  absolute  value  of  the  range  error  after  correction,  i.e., 
the  difference  between  the  moduli  of  zc  and  z.  Finally,  since  the  argument 
of  zc  is  the  same  as  that  of  zc,  the  bounds  of  figure  10  apply  to  the 
absolute  value  of  the  angle  error  after  correction,  i.e.,  the  difference 
between  the  arguments  of  zc  and  z. 

From  the  standpoint  of  implementation,  correction  might  best  be  accomplished  by 
means  of  a polynomial  approximation  to  (48).  Along  the^e  lines,  the  Taylor 
Remainder  Formula  can  be  used  to  show  that 
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i 


p 


[1  + H/2E  - p2/8E2  + e (p,  H,  E)] 


(52) 


f(P,  H,  E) 


P ~ 


\Jl  + K/E 


where 


I e (P,  H,  E)|  < 1.5  x 10-6 

whenever  p < 200  nautical  miles,  H < 60,000  feet,  and  E satisfies 
(49).  Thus  letting 

f (p,  H,  E)  - p - p [1  ♦ H/2E  - P2/8E2] 

o ______ 

Vl  + K/E 


(53) 


(54) 


error  correction  can  be  accomplished  by  adding  fQ  (Pc,  H,  E)  to  pc  without 
a significant  change  in  the  final  outcome  due  to  approximation  error. 

As  pointed  out  in  Section  6,  zc  is  a function  of  target  range  P,  target  angle 
Y , and  target  altitude  H.  Consequently,  the  conversion  error  after  correction 
with  the  cubic  approximation  (54)  can  be  viewed  as  a function  of  these  three 
parameters  and  E,  i.e., 

6[  (p,  Y,  H,E)  -|  « - (pc  ♦ fD(pc,  H,  E)]e  lYc  I . (55) 

The  maximum  of  this  error  with  respect  to  Y is  illustrated  in  figure  17  for  a 
radar  at  sea  level  on  parallel  38  and  a target  at  60,000  feet.  The  value  of 
the  parameter  E has  been  taken  to  be  the  polar  radius  of  the  reference  ellip- 
soid in  this  example.  Actually,  the  conversion  errors  after  correction  shown 
in  figure  17  are  greater  than  those  observed  for  targets  at  altitudes  less 
than  60,000  feet.  Moreover,  as  expected,  in  all  cases  the  errors  are  lower 
bounded  by  the  range  error  oscillation  amplitude.  Conversion  errors  after 
correction  corresponding  to  the  case  where  the  radar  site  altitude  is  10,000 
feet  are  essentially  identical  to  those  obtained  in  the  situation  where  the 
radar  is  at  sea  level.  For  this  reason,  the  altitude  of  the  radar  site  does 
not  appear  in  figure  17.  Hence,  figure  17  can  be  viewed  as  representing  upper 
and  lower  bounds  on  conversion  error  after  correction  for  tanget  altitudes  not 
exceeding  60,000  feet.  Insofar  as  the  size  of  conversion  error  is  concerned, 
a comparison  of  figure  17  with  figures  14  and  15  indicates  that  error  correction 
offers  some  advantages  as  an  error  control  device  that  cannot  be  realized  with 
the  minimax  method. 

While  the  correction  approach  to  error  control  does  seem  to  offer  some  advan- 
tage over  the  minimax  method,  the  practicality  of  the  procedure  is  still  a moot 
question.  For  example,  we  have  compared  the  two  methods  under  the  assumption 
that  slant  range,  altitude,  and  azimuth  are  known  exactly.  It  remains  to  consider 
the  situation  in  which  there  is  measurement  error.  In  addition,  we  have  restricted 
ourselves  to  the  case  where  K is  61,000  feet.  In  the  face  of  measurement  error, 
this  may  not  be  the  best  choice  for  K.  Finally,  there  may  be  other  limitations 
imposed  by  hardware  and  software  considerations. 
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FIGURE  17.  MAXIMUM  CONVERSION  ERROR  AFTER  CORRECTION  FOR  RADAR  ON  PARALLEL  38 
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CONCLUDING  REMARKS 
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Numerical  results  have  been  obtained  for  conversion  errors  in  the  context  of 
an  earth  model  based  on  the  reference  ellipsoid.  Radar  locations  between 
parallels  28  and  48  at  altitudes  between  0 and  10,000  feet  were  considered. 
Attention  was  directed  to  targets  at  altitudes  from  1,000  to  60,000  feet. 
Differences  between  geodetic  and  conformal  latitudes  were  taken  into  account. 
Results  have  been  presented  in  terms  of  ranges  and  angles  measured  in  the 
local  stereographic  projection  plane  with  the  projection  of  the  radar  site 
location  taken  as  the  origin  of  coordinates.  Target  ranges  from  5 to  200 
nautical  miles  were  considered. 

At  constant  target  range  and  altitude,  both  conversion  range  and  angle  errors 
can  be  viewed  as  oscillatory  functions  of  target  angle.  The  median  of  an 
angle  error  oscillation  is  0,  whereas  the  median  of  a range  error  oscillation 
is  dependent  upon  target  range  and  target  altitude.  The  amplitude  of  an  angle 
error  oscillation  is  invariant  to  radar  site  altitude  and  the  local  projection 
plane  height.  This  is  also  essentially  true  in  the  case  of  the  range  error 
amplitude.  In  addition,  its  dependence  on  target  altitude  is  negligible.  On 
the  other  hand,  both  amplitudes  increase  monotonically  with  target  range,  and 
angle  error  amplitude  is  essentially  a nondecreasing  function  of  target 
altitude.  For  radar  site  locations  between  parallels  28  and  48  at  altitudes 
from  0 to  10,000  feet,  results  indicate  that  angle  error  amplitude  is  upper 
bounded  by  a linear  function  of  target  range  assuming  the  values  of  .00025* 
and  .0056*  at  target  ranges  of  5 and  200  nautical  miles,  respectively.  By  the 
same  token,  range  error  amplitude  is  upper  bounded  by  a convex  up  curve  that 
is  essentially  0 at  a target  range  of  5 nautical  miles  and  increases  to  0.0195 
nautical  miles  at  a target  range  of  200  nautical  miles. 

The  height  of  the  local  stereographic  projection  plane  was  determined  by  two 
methods.  One  of  these  corresponds  to  a minimax  approach  to  error  control 
based  on  a spherical  earth  model.  The  other  gives  rise  to  another  approach 
for  reducing  conversion  error.  In  the  case  of  the  minimax  procedure,  the 
median  of  a range  error  oscillation  is  invariant  to  radar  site  latitude.  On 
the  other  hand,  it  is  strongly  dependent  upon  the  altitude  of  the  site.  In 
the  case  of  the  alternative  approach  to  error  control,  it  is  essentially 
invariant  to  all  aspects  of  radar  site  location  thereby  providing  a set  of 
universal  error  correction  curves. 

In  the  case  of  the  minimax  procedure,  the  overall  conversion  error  can  reach 
0.08  nautical  mile  for  a radar  at  sea  level  on  parallel  38.  If  the  site  is 
raised  to  10,000  feet,  errors  as  large  as  0.22  nautical  mile  are  possible. 

In  the  case  of  the  alternative  approach  to  error  control,  results  indicate  that 
the  overall  conversion  error  can  be  reduced  considerably.  In  fact,  numerical 
results  indicate  that  the  overall  conversion  error  is  essentially  reduced  to 
the  level  of  the  range  error  oscillation  amplitude.  Cubic  approximations  to 
error  control  curves  derived  from  this  approach  appear  to  be  capable  of 
restricting  the  overall  conversion  error  to  within  2 meters  of  this  level. 
However,  the  practicality  of  the  approach  is  a moot  question. 
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In  conclusion,  it  must  be  emphasized  that  we  have  considered  conversion  errors 
under  the  assumption  that  sensor  measurements  are  exact  and  data  processing  is 
carried  out  with  infinite  precision.  In  practice,  this  is  not  the  case. 
Receiver  front-end  noise  exists,  range  and  azimuth  measurements  are  quantized 
before  being  passed  on  to  the  central  processor,  trigonometic  functions  and 
other  mathematical  relationships  are  approximated  in  various  ways,  and  multi- 
plication, addition,  etc.,  are  carried  out  using  finite  precision  arithmetic. 
All  of  these  contribute  to  error,  and  we  have  no  reason  to  believe  that  the 
effect  upon  the  conversion  process  is  inconsequential.  In  this  sense,  our 
reiulti  might  be  viewed  as  representing  lower  bounds  for  conversion  error  that 
•re  attainable  only  under  ideal  operating  conditions. 
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